function u = func_ex1(n, m, i, j)
    global bounds
    
    x = 0 : bounds(1) / (n - 1) : bounds(1);
    y = 0 : bounds(2) / (m - 1) : bounds(2);
    
    u = x(i)^2 + y(j)^2;
end

